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Abstract: While chiral symmetry breaking in the QCD vacnnm is attribnted to nonzero 
chiral condensate, an alternative symmetry breaking pattern with no chiral condensate is 
also possible, as pointed ont by Stern. This hypothetical phase was exclnded in QCD at 
zero density long time ago, bnt nothing forbids it at hnite baryon density. In this work, we 
stndy the 9 dependence of this nnorthodox phase on the basis of chiral pertnrbation theory. 
Physical observables snch as energy density, topological snsceptibility, non-local chiral order 
parameter and meson masses are compnted analytically in the e-regime. At nonzero 6 we 
hnd an exotic phase that breaks vectorial flavor symmetries in a way analogons to the Aoki 
phase in lattice QCD. 
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1 Introduction 

Spontaneous symmetry breaking can be characterized by order parameters that transform 
nontrivially under the symmetry of interest. A commonly used order parameter for chiral 
symmetry breaking in Quantum Chromodynamics (QCD) is the chiral condensate, 

It is linked to the accumulation of near-zero Dirac eigenvalues through the Banks-Casher 
relation [ 1 ]. Another order parameter is the pion decay constant, F-jr. Some time ago 
Stern [2, 3] pointed out that the condition for 7 ^ 0 is weaker than that for (V’V’) 7 ^ 0) 
suggesting the possibility of an exotic phase in QCD in which = 0 but 7 I 0.^ We 

will refer to this phase as the Stern phase. Chiral symmetry breaking in this phase could 
be triggered by four-quark condensates such as + 75)'*/’) and 

^?/)A“(l — 75 )';/’• + 75 )^^)) with {A“} the flavor generators [4].^ These condensates 

leave the discrete anomaly-free subgroup of U(l)y 4 unbroken, which ensures a vanishing 
chiral condensate. (Actually the possibility of an unbroken discrete axial symmetry was 
pointed out by Dashen long time ago [6].) The Stern phase is analogous to antiferromagnets 
which has no global magnetization; also similar is the so-called molecular Bose-Einstein 
condensation (BEC) [7, 8] which is distinguished from the atomic BEC by an unbroken Z 2 
symmetry. 

^Other bilinear condensates such as are assumed to vanish as well. 

^While quartic condensates also form in color-superconducting phases of QCD at high density [5], the 
baryon number symmetry U(l )b is not broken in the Stern phase. Moreover the patterns of U(1 )a symmetry 
breaking in the Stern phase and in the color-superconducting phases are different, as will be explained later. 
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Soon after the proposal, the Stern phase in QCD was critically examined in [4] where 
it was proved with rigorous QCD inequalities that this phase is ruled out in QCD at any 
temperature and zero density. This proof, however, leaves open the possibility that the 
Stern phase may emerge in QCD at nonzero chemical potential since the complex path- 
integral measure invalidates the use of QCD inequalities. Indeed, a Ginzburg-Landau-type 
analysis in a chiral effective model suggests that this is likely to be the case [9]. In addition, 
studies of inhomogeneous chirally broken phases in dense QCD suggest that chiral symmetry 
breaking in such phases could be driven not by the chiral condensate but rather by a higher- 
order condensate [10, 11]: the basic idea of [11] is that a one-dimensionally modulated chiral 
condensate is wiped out by thermal fluctuations of phonons, whereas [10] shows in QCD 
at large N that a higher-order chiral order parameter whose spatial average is nonzero 
must exist when the chiral condensate is locally nonzero but its spatial average vanishes. 
Recently, phases with massive fermions with no bilinear condensate have been found in 
numerical simulations [12, 13], which bears resemblance to the Stern phase. So, even 
though the presence of the Stern phase in QCD remains an open problem for now, we have 
pieces of circumstantial evidence suggesting that the Stern phase is a realistic possibility 
worthy of serious consideration. This will be of importance for our better understanding of 
the QCD phase diagram, which is still only poorly understood [5, 14]. 

If the Stern phase indeed exists in the finite-density QCD, there must be a transition 
from a hadronic phase to the Stern phase as /i is varied. A possible phase structure at 
// 7 ^ 0 was proposed in [9] where two transitions were reported along the ^ axis: from 
the hadronic phase to the Stern phase and then to the chirally symmetric phase. These 
can become smooth crossovers for nonzero quark masses. At the first transition the chiral 
condensate drops dramatically while the four-quark condensate is unaffected. A rapid rise 
of the baryon number susceptibility is a signal of this transition. We would also like to 
mention another scenario based on the idea of inhomogeneous condensation [11] in which 
the above transitions are both second order in the chiral limit; at the low-/x transition, it is a 
proliferation of domain walls and associated Nambu-Goldstone modes that drive the chiral 
condensate to zero. These pictures are based on effective models and a quantitative precision 
is not expected, but nonetheless their symmetry-based arguments are robust predictions 
that can be tested in QCD-based calculations in future. 

In this paper, we investigate various aspects of low-energy physics in the Stern phase by 
means of chiral perturbation theory. In particular the structure of the 9 vacua in the Stern 
phase is analyzed in great details for the first time. We find behaviors that differ drastically 
from those in the orthodox 9 vacuum. The competition between multiple leading terms in 
the chiral Lagrangian is shown to lead to a nontrivial phase diagram at nonzero 9. Not 
only quarks in the fundamental representation of the gauge group, but also those in higher 
representations are considered and new results are obtained. 

This paper is organized as follows. In section 2 we sort out the breaking pattern of 
continuous and discrete symmetries in the Stern phase and present a systematic derivation 
of the chiral effective theory. While this part overlaps with preceding works [4, 15, 16], 
we extend them by considering the most general breaking pattern of the discrete axial 
symmetry. One of the new results here is the existence of topologically stable domain walls 
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in the Stern phase, and another is an analytical calculation of the volume dependence of a 
non-local chiral order parameter in the e-regime. In section 3 we introduce the 0 angle into 
the low-energy effective theory and compute various observables such as energy density, 
topological susceptibility, topological density and pion masses. It is revealed that low- 
energy physics at 0 7 ^ 0 is sensitive to a subtle balance between leading terms in the chiral 
Lagrangian. An exotic phase similar to the Aoki phase of Wilson fermions [17] is uncovered, 
and its domain of existence is determined in the phase diagram. Finally we consider QCD 
with quarks in higher representations and elucidate a multi-branched 6 dependence of the 
energy density that surprisingly differs from the case of fundamental quarks. We conclude 
in section 4. 

2 Chiral effective theory 

2.1 Symmetries and effective Lagrangian 

In this section we classify low-energy chiral effective theories for the Stern phase, general¬ 
izing preceding works in [4, 9, 15, 16]. We will also give a brief account of topologically 
stable domain walls in the Stern phase, which has not been discussed to date. 

2.1.1 Massless quarks 

Let us consider SU(A^) gauge theory with Nf > 2 massless Dirac fermions in a complex 
representation R of SU(A^) in Euclidean spacetime. As is well known, the classical U(1 )a 
symmetry in the chiral limit is violated by quantum effects due to instantons [18, 19], but 
generally there exists a discrete remnant of the U(1 )a symmetry. According to the index 
theorem, the index Ir of the Dirac operator in the representation R for a single instanton 
background is given by [20, 21] 

Ir = ‘^Tr , (2.1) 

with Tr defined by tr(r“T^) = Tr5°'^ for SU(A^) generators in the representation R\ e.g., 

1 N + 2 N -2 

Tpund = 2 > ^Adj = Ai, Ts = —-— and Tas = —-— , (2 .2) 

where S (AS) stands for the two-index symmetric (anti-symmetric) representation of S\J{N), 
respectively. This implies that the U(1 )a symmetry shrinks to 'Z^Nj.Tji due to quantum 
effects. Then the orthodox pattern of chiral symmetry breaking with 7^ 0 reads^ 

S\]{Nf)R X SU(A^/)l X U(1)b X {Z^N,Tn)A S\]{Nf)v x U(l)s • (2.3) 

By contrast, the putative Stern phase entails a different pattern of chiral symmetry breaking: 

S\]{Nf)R X S\]{Nf)L X U(1)b X {Z^NfTn)A S\]{Nf)v X U(1)b x {Zk)a , (2.4) 

^Precisely speaking, Z2jvy C J-uNjTfi is part of the SU(W)^? ^ SU(W)i ^ U(1)b group and has to be 
factored out in the LHS of (2.3) to avoid double counting. For the same reason, Zjvy. which is contained 
in both S\J{Nf)v and U(l)s has to be factored out in the RHS of (2.3). However we shall be cavalier on 
these formalities for simplicity of exposition. 
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where 2 < K < ANjTr is an even divisor of (Note that, when K = 2, (2.4) 

is equivalent to (2.3) and there is nothing new.) While only the specihc case of K = 
dA^jTpund = has been discussed in the literature [4, 9, 16], the other values of K are 
also theoretically admissible. As the residual {J-k)a with K > 2 enforces = 0, chiral 

symmetry breaking in the Stern phase must be driven by higher-dimensional condensates 
[4]. For example, {{'tpRipi)"') + h.c. 7 ^ 0 for some n>2 corresponds to iF = 2n, whereas 
(V'ijV’L • 7 ^ 0 or ([det('i/’^V'£)]^'^^) 7 ^ 0 (this is nothing but the ’t Hooft vertex) 

corresponds to the maximal unbroken symmetry, K = AN 

The vacuum structure of the Stern phase must be understood with some care. Since 
(I-ANfTji)A is spontaneously broken to {J.k)a, it appears at hrst sight that there will be 
ANfTji/K isolated degenerate vacua. This is not quite correct, however. The point is that 
two vacua that can be rotated to each other via an action of (I- 2 Nf)A are not isolated, but 
are continuously connected to each other with no potential barrier via a non-Abelian chiral 
transformation. This is obvious from the fact that (Z 7 Vy^)_R x (ZAr^)^ C S\J{Nf)jixS\J{Nf)L- 
As a result, the would-be domain walls separating such vacua are unstable, as stressed in 
[22] for a fractional axial domain wall in the QCD vacuum. 

Then, under what conditions does a stable domain wall exist in the Stern phase? 
Evidently there must be multiple vacua that cannot be rotated to each other via a combined 
action of (Z^c)^ and (I- 2 Nf)A- It is not difficult to see that this is true if and only if 

LCM(A:, 2Nf) < ANfTn , (2.5) 

where LCM(a, b) for a, 5 G IM is the least common multiple of a and b. We also see that 

( the number of disconnected components 
of the vacuum manifold 

The bottom line is that a stable domain wall can exist under the condition (2.5) and that 
the variety of domain walls is determined by (2.6). Let us make a few quick comments. 
First, the RHS of (2.6) is always a positive integer because both K and 2Nf are divisors 
of ANfTji. Secondly, (2.5) cannot be satished if Tr = 1/2, implying that the vacuum 
manifold is connected for fundamental quarks for any K. Thirdly, when K = 2 (i.e., the 
QCD vacuum with chiral condensate), the number of isolated vacua is given by 2Tji, as 
follows from ( 2 . 6 ). 

Let us comment on the literature. It is well known that in A? = 1 SU(A^) Super- 
Yang-Mills (SYM) theory, 'Z 2 N C U(l)ij breaks down spontaneously to Z 2 through gaugino 
condensation [23, 24]. There are N isolated ground states that are discriminated by phases 
of the condensate as {XX) ~ A^ exp(27rzfc/A^) with A; = 0,1,..., — 1, and stable domain 

walls exist [25-27]. We also wish to mention the so-called axion domain walls [28], which 
have been discussed widely in axion cosmology. There is an apparent similarity between 
domain walls in these theories and those in the Stern phase, and we anticipate that many 
properties would be shared in common. Nevertheless, it deserves attention that the Stern 

should be even because (Z 2 )a that flips the sign of ip is just a 27r rotation of space, which cannot be 
spontaneously broken when the Lorentz symmetry is intact. 


^NfTR 


LCM{K, 2Nf) 


( 2 . 6 ) 
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phase possesses gapless pion excitations, which are missing in SYM and the axion theory. 
It would be intriguing to explore physical consequences of this difference in details. 


2.1.2 Massive quarks 

Next we switch on the mass term iii th® microscopic Lagrangian, with M 

the NfxNf quark mass matrix. This term breaks {J.k)a down to (I- 2 )a explicitly. However, 
one can make the Lagrangian invariant under (J.k)a if M along with ipn/i transform as 

and M ^ M for ^^^^{Tk)a- (2.7) 

This symmetry should be preserved in the low-energy effective Lagrangian £(?7, M) of the 
N‘j — 1 Nambu-Goldstone modes, pions, denoted collectively by U{x). Since the axial 
current is neutral under {'I.k)a, so are pions. Thus the condition on the 

effective theory imposed by {'I.k)a invariance reads 

£([/, M) = C{U, zM) for E Zk /2 ■ (2.8) 

Using the Nf x N^ coset variable U E SU(A^j)a as a building block, one can straightfor¬ 
wardly write down the most general chiral Lagrangian consistent with (2.8) [16]. Let us 
start the classification with K = A and Nf > 2. Then all the odd powers of M are forbidden 
by (2.8) and we find, up to second order in d and M, 

^NT> 2 iU,M) = {OaU^ d^U + v^diU^diU) - /i|tr(M[/)| 

{tv{MU)f + tr{MU)^ ^ 

- gi tr(MMi). (2.9) 

The low-energy constants {h,hi,h 2 ) and the high-energy constant gi are analogous to Lg, 
L 7 , Lg and H 2 in standard chiral perturbation theory at 0{p^) [29]. They are related to 
chiral susceptibilities, hi and h 2 can be complex in general. The pion velocity v can differ 
from the speed of light owing to the breaking of Lorentz symmetry in medium. We note 
that /, u, h, hi, /i 2 and gi all depend on the chemical potential g implicitly. 

For two flavors, the identity — tr(M?7)^ = 2detM allows us to cast the 

leading-order Lagrangian in the form 

f2 

£^=4(77,M) = ^tr {d^U^d^U + v^diU^diU) - h|tr(Mt/)| - {hti:{MUf + h.c.} 

— g'l tr(MM'*') — (52 det M-|-h.c.). (2.10) 

While the last term is independent of the pion field, it is 0-dependent and contributes to 
the topological susceptibility (cf. section 3.1). As a side remark we mention that it plays 
an important role in QCD at high temperature [30]. 

The absence of the linear term tr(M[/) -|- h.c. in (2.9) and (2.10) is consistent with 
= 0 in the chiral limit. This implies that the power counting in the p-regime of this 




(tr(MU)) - tr(M[/)2 


-I- h.c. 
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phase must be modified from the usual one, d r\j 0{p) and M ~ 0{p^), to 9 ~ M r-u 0{p).^ 
As for the pion mass, (2.9) and (2.10) imply tUtt oc M [4], in contradistinction to the 
conventional picture where rn^ oc M. Actually the abnormal scaling ttItt oc M has been 
known for the color-flavor-locked (CFL) phase of dense QCD [32-34], superfluid phase of 
dense two-color QCD [35] and two-flavor QCD at 9 = tt [36-38]. What is common in all 
these cases is that the chiral condensate is either zero or negligibly small. 

The three 0{M'^) contributions in (2.9) have different origins. If the order parameter 
of symmetry breaking is in the adjoint representation of the leading mass- 

dependent contribution should be trAdj(AfC/) = |tr(MC/)|^ — 1. Thus the term oc /i in (2.9) 
originates from a condensate that transforms in the adjoint flavor representation. Similarly, 
the term oc hi (oc / 12 ) comes from a condensate in the two-index symmetric (anti-symmetric) 
flavor representation, respectively.® 

Let ns finally consider K > 47 This time the effective theory is considerably simplified: 

f2 

= ^(d4l7'^d4U+ - /i|tr(M[/)| - c/itr(MMt) (2.12) 

at leading order. The other terms are banned by the discrete symmetry (2.8). 

Consequently, enjoys invariance under an arbitrary phase rotation of M, which is 

equivalent to a vanishing topological susceptibility at this order. In other words, topologi¬ 
cally nontrivial sectors are entirely suppressed. 

The classification of the effective theory [(2.9), (2.10) and (2.12)] for general K is the 
main result of this subsection. 

A brief comment is in order concerning hadrons other than pions. For the effective 
theory of pions to be a valid low-energy description, the baryon sector must have a mass 
gap.® However the status of baryons in the Stern phase is still elusive. In principle, 
dynamical masses of fermions can be generated without bilinear condensate as evidenced 
in [12, 13]. Also in QCD, it has been recognized historically that baryons in the “mirror 
assignment” can acquire a dynamical mass even when the chiral condensate vanishes [42- 
44]. Previous researches on the Stern phase have found that the baryon spectrum in the 
Stern phase crucially depends on the chirality assignment of baryons [4, 9]. Instead of 
trying to resolve this delicate issue, we shall content ourselves in this paper by assuming 

®This is similar to the generalized chiral perturbation theory [31]. Note however that Ref. [31] retains 
nonzero chiral condensate whereas we put it to zero exactly as a consequence of the (Zk)a symmetry. 

®In the CFL phase, the leading gauge-invariant order parameter of chiral symmetry breaking is the 
four-quark condensate [39, 40] 

{efgheijkipRtp^ifililii) oc Shk , (2.11) 

which transforms in the anti-symmetric representation of SU(3)i{/L. Reflecting this pattern, the chiral 
Lagranglan in the CFL phase only contains the term oc /12 in (2.9); h = hi = 0 is indeed conhrmed in 
explicit microscopic calculations [33] . 

^We note that the condition if > 4 is consistent with K < ANfTu only if NfTji > 1. This means that 
R must be higher than fundamental for Nf = 2, whereas no such constraint arises for Nf > 2. 

®Since chiral symmetry is spontaneously broken in the Stern phase, the ’t Hooft anomaly matching 
condition [41] is satished by pions. If some baryons happen to be gapless, then it would be highly nontrivial 
to keep the anomaly matching satisfied. This seems to be a rather unlikely possibility. 
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that the baryon sector is fully gapped. It should also be kept in mind that, at high density 
or in the large-iV limit, instanton effects are suppressed and the pseudo-Nambu-Goldstone 
mode {r]') associated with U(1 )a breaking gets light. Then one has to incorporate rj' into 
the effective Lagrangian as well. This is an intriguing situation but will not be covered in 
this paper. 


2.2 Finite-volume partition function 

When pions are sufficiently light in a finite volume, their zero-mode fluctuations become 
non-perturbative and have to be integrated out exactly. This occurs in the so-called e- 
regime [45, 46] where the linear extent L of the Euclidean box is such that 

—. (2.13) 

AqcD 'THtt 

This means that the contribution of hadrons other than pions to the partition function 
is suppressed (first inequality) whereas pions’ Compton length is sufficiently larger than 
the box size (second inequality) so that the non-zero modes of pions become irrelevant. 
This regime can be realized by taking the double limits V 4 = —)• oo and M —)• 0 with 

fixed. More formally stated, we shall adopt an exotic e-expansion scheme with 
d 1 '^ l/L 0(e) and M ~ O(e^) [15]. A similar scheme was used in dense QCD [35, 47] 
but it differs from the conventional e-expansion with M ~ 0{e'^). This disparity of course 
stems from the absence of the linear mass term in the Stern phase. 

In this limiting regime, the QCD path integral reduces to a finite-dimensional integral 
over pion zero modes. For the three cases in section 2.1 we obtain, respectively,® 


yK=4: 

^Nf>2 


(M) 


z§;Mm) 


'SV{Nf) 


dU exp 


K. ( hVr(Mut + (+1,,; 








(2.14a) 


dU exp F4{/i|tr(M17)|^+ (htr(M[/)2 + h.c.) +giiv{MM^) 


'SU(2) 

-|- (52 det M + h.c 


>} 




dC/ exp F 4 {/i|tr(MC/)p + 5 itr(MMt)| 


(2.14b) 

(2.14c) 


where dU denotes the Haar measure. These expressions give exact mass and volume- 
dependence of the partition function in the ^-regime, which is not only theoretically inter¬ 
esting but also useful in that we can extract physical quantities in the infinite-volume limit 
(e.g., low energy constants h,hi and / 12 ) from numerical data obtained in a finite volume 

®A similar finite-volume analysis was performed in [15], but the authors did not specify the discrete 
symmetry responsible for the vanishing chiral condensate, nor did they underline the distinction between 
AT = 4 and K > 4: whose importance is clear from (2.14). Quarks in higher representation of the gauge 
group were not considered, either. 
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through fitting to finite-volume formulas extracted from (2.14).^^ We hope that analytical 
results in this section serve as a guide in future lattice simulations of the Stern phase. 


• Partition function for K > 4 Since the structure of the partition functions (2.14) 
is mathematically more involved than in the conventional e-regime, we shall focus our 
attention on the K > A case, (2.14c), for simplicity. As S\J{Nf)v symmetry is assumed to 
be unbroken in the Stern phase, we assume that U = 1 is the ground state. This fixes the 
sign of h to positive. Extending the manifold of integration to \J{Nf), we obtain 


exp 


-P45itr(MMt)j 


lv{Nf) 


dU exp 


= - I dh e-l"l' 

TT Jc 

roo 

= 21 dx xe-^" 

'o 


V 4 htr(MU) tr([/tMt) 

dU exp \z^/vJiiT{MU) + z*M^) 




\ r dx xe-^' / 
Jn J\] 


dU exp 


V{Nf) 


2x^/vji Retr(M[/) 


(2.15) 

(2.16) 
(2.17) 


In the last step the phase of z was absorbed in U. We now set M = diag {ruf) and dehne 
/xj = 2y/V^ ruf. Assuming G IR we substitute the well-known analytic formula for the 
above unitary integral [49-51] to obtain 


Z 


K>4: 


({/^/}) = CNf exp 


Ah 




dxxe~^" , (2.18) 


f / A((x/xi)2,...,(x/xjvp2) 

where ..., aj^) = ]^(®?~Oj) is the Vandermonde determinant, and the normalization 

W /-1 


i>j 


constant C^f = i)/2+i 


[ n 


k=l 


ensures —)• 1 in the chiral limit. In particular. 


for Nf = 2 and /xi = /X 2 = /x we have 

= exp 


91 +h 2 








(2.19) 


Equations (2.18) and (2.19) are new results. Let us contrast (2.19) with the = 2 
partition function in the topologically trivial sector of the conventional e-regime [46]: 

Zjy=o({w., m}) = Iq (V451m) — if (I^Sm). (2.20) 

• Chiral susceptibility Although chiral condensate in the Stern phase vanishes in the 
chiral limit, there is a non-local order parameter for chiral symmetry breaking. Considering 
Nf = 2 for simplicity, we define the disconnected chiral susceptibility 


Xudim) = / d^x {uLUR{x)dRdLiO)) - {urur) (dRdi) 

1 52 


-|- h.c. 


= lim 


V4 druudm'* 


log Z -|- h.c., 


( 2 . 21 ) 

( 2 . 22 ) 


refer the interested reader to [48] for a review of other miscellaneous intriguing aspects of QCD in 
a finite volume. 
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Figure 1. The non-local order parameter (2.21) in the microscopic limit for itT > 4 and Nf = 2 
with equal masses. 


which is singlet under a vectorial isospin rotation but is charged under the axial isospin 
rotation generated by 75 ra. Thus Xud 7 ^ 0 in the chiral limit is a signal of spontaneous 
breaking of SU(2)^. Noting that giti{MM^) = 91 no contribution 

to Xud, we have, for it' > 4 in the e-regime 


1 52 If 

Xud{m)= lim — --log / dU exp 

mu,d^rn V4 dniuOm*^ \JsV{2) 


1 / 4/1 tr(Mi/) tr(i/^Mt) 


+ h.c. 


(2.23) 


This expression can be evaluated analytically. After a tedious calculation one finds 

,2 / r r N 21 


Xud — h 


2 Ii — I 2 ^ 9? Iq + h — h — h 


3 io — ii 


6 


/ 0-/1 


9 ^ ( h-h 
h-h 


(2.24) 


where in’s are the modified Bessel functions of the first kind, In{x), evaluated at x = 9 ^I2. 
The derivation of (2.24) is lengthy and is relegated to appendix A. Figure 1 shows Xud as 
a function of /r^/2 = 2 V 4 hm?. Asymptotically Xud behaves as 

~ ^9^ for /i <C 1 and ~ 1 -|-^ for /r » 1. (2.25) 

The vanishing of Xud for /r —)• 0 is indicative of symmetry restoration at finite volume. In 
the opposite limit /i 1 , Xud approaches /i, so /i 7 ^ 0 is in fact an order parameter of chiral 
symmetry breaking, much like S in the conventional e-regime. It is quite intriguing that 
the behavior of Xud in figure 1 is non-monotonic: it approaches h from above! This feature 
is not seen in the volume dependence of chiral condensate in the QCD vacuum [45]. 


• Spectral sum rules The mass dependence of the partition function also provides 
detailed information on the statistical distribution of Dirac eigenvalues. Let us first observe 
that the QCD partition function in the topologically trivial sector may be cast in the form 


■^QCD 
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(2.26) 













where {i\n}n denotes eigenvalues of the Euclidean Dirac operator and the primed product 
runs over eigenvalues with ReA„ > 0. The bracket represents the expectation value 

of O with the weight of QCD with Nj massless flavors. By expanding .^qcd in 
and equating the coefficients with those from the effective theories (2.14a), (2.14b) and 
(2.14c), one obtains an inhnitely many spectral sum rules obeyed by Dirac eigenvalues.^^ 
When K > A, the spectral sums generally depend on both gi and h. However, the term 
gi tr(MMl) appears even in a free theory and has no bearing on symmetry breaking at low 
energy. Rather, it serves to absorb UV divergences arising from large perturbative Dirac 
eigenvalues [15, 52]. Therefore we should make suitable combinations of spectral sum rules 
in such a way that gi does not appear explicitly. In terms of rescaled dimensionless Dirac 
eigenvalues Cn = 2^/14^ A^, the first few sum rules obtained this way for K > A read 




(2.27a) 

(2.27b) 


The existence of such nontrivial correlations on the scale Cn ~ ^>(1) suggests that the 
typical scale of Dirac eigenvalues pertinent to symmetry breaking in the Stern phase is 
~ l/y/V^h. This volume dependence is exactly in accord with the prediction by Stern 
[2, 3]. While this is in contrast to the conventional microscopic domain of the QCD vacuum 
where An ~ I/V 4 S, there is a similarity to the microscopic domain in high-density QCD 
where An ~ 1/\A44^ with A the BCS gap of quarks near the Fermi surface [53], indicating 
a natural correspondence -H- h. Finally we point out that the reality of the spectral 
sums (2.27) is rather nontrivial, because the Dirac operator is not assumed to be anti- 
Hermitian and eigenvalues (^n are complex-valued in general. This may be pointing to a 
hidden symmetry in the Dirac spectra of the Stern phase. 


3 The 6 vacua 

Physics of the 9 vacuum in QCD has been investigated over many years. Not only is it 
relevant to the so-called strong CP problem and axion physics, it has recently gained a 
renewed interest in the context of possible CP violation in heavy ion collisions [54-56]. On 
a practical side, lattice simulations at fixed topology suffer from large hnite-volume effects 
and it is useful to analytically understand the topology dependence of observables [57, 58]. 
Various aspects of ^-dependent physics are reviewed in [59]. 

While 0-dependence in QCD is inherently nonperturbative, it is quite difficult to simu¬ 
late QCD with a nonzero 9 angle on the lattice because of a severe sign problem. It is then a 
promising alternative to employ chiral effective theories to study topological aspects of QCD 
at low energy [60-63] . Studies of the 9 vacuum in the e-regime of chiral perturbation theory 

^^When Tf > 4 all the topologically nontrivial sectors are gone, while for K — 4 the mass expansion in 
each topological sector leads to different sum rules. 
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were performed by various authors [46, 64-66]. The virtue of taking the microscopic limit is 
that the partition function and other various nonperturbative quantities (e.g., topological 
susceptibility) can be computed exactly. In this section we shall extend this analysis to the 
Stern phase and delineate the structure of the 6 vacua, emphasizing qualitative differences 
from the ordinary QCD vacuum. 


3.1 Quarks in the fundamental representation 

The partition function in the presence of the 0 angle admits a Fourier decomposition 

OO 

z{e)= ( 3 - 1 ) 

Q=—oo 

where Zq is the partition function in the sector of topological charge Q with 

Q = ^ / • (3.2) 

When there are Nj quarks in the fundamental representation, the 9 angle can be transferred 
to the complex quark mass matrix as M —>■ M via an axial rotation. As is evident from 

(2.14), Z^^‘^{M) has no dependence on the 9 angle: at leading order in the e expansion, 
Zq with Q 7 ^ 0 do not contribute and we simply have Z{9) = Zq. 

We now focus on the K = A case, and especially Nf = 2 for simplicity. 

• Partition function with the 0 angle Substituting M = in (2.14b) yields 


Zf=i2(m,0) = 


/SU(2) 


dU exp 


V 4 m^ |/i(tr 17)^ + 2gi + (he^^ tr(17^) + h.c.) + [§2 e*^ +h.c.)| 

(3.3) 


Although 9 dependence is strongly affected by the phases of h and g 2 we currently lack 
information on their physically appropriate values. To get an idea of how Z depends 
qualitatively on 0, let us assume that h and g 2 are both real. Then, using the identity 
(tr U)‘^ — tr([/^) = 2, one obtains 


^Nf =2 


(m, 9) 


_ ^2V4m'^[gi + {g2-2h)cose] 



^V4,m?{h+2h cos 6 )(tr U)^ 


^ ^2V4mHgi+h+g2Cos9) - Ii{ae)} 


a.0=2V4 {h-^2h cos 


(3.4) 

(3.5) 


where in the last step we have used (A.6a) in appendix A. This is the starting point of our 
analysis in this subsection. We demand that U = 1 he the ground state at 9 = 0, which 
translates into the condition 


h + 2h>0. (3.6) 

To illustrate the physical content of (3.5) we shall calculate two quantities of major physical 
interest: topological susceptibility xt and the energy density E{9). 
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3.5 



Figure 2. Convergence of the topological susceptibility Xt in the microscopic limit toward Xt° 
K = 4 and Nf = 2 with equal masses. 


• Topological susceptibility Let us recall the definition 


Plugging (3.5) into this dehnition, we obtain 


0=0 


M 

^4 


Xt = ‘^rn^ \g 2 -h 


loi^o) — 2 / i ( q ! o ) + -^ 2 ( 00 ) 1 


-^o(«o) — hioto) 


7 


ao=2V4(h-\-2h)m‘^ 


(3.7) 


(3.8) 


In particular, in the macroscopic limit (oq 1 )/^ one hnds xt xV 2 ut ,^(52 + 2/i). 
While it is natural that Xt^ vanishes in the chiral limit, we hnd it interesting that xT ^ 
in contrast to the conventional behavior xt = '^rn/Nf oc m in the QCD vacuum [46]. Figure 
2 plots the deviation of Xt from as a function of the scaling parameter oq- We note 
that it is the combination h + 2h that controls the hnite-volume effect for xt and that gi 
and g 2 play no role here because they do not couple to the pion fluctuations. 


• Energy density Next we calculate the 0-dependent energy density dehned by 

E{e) = -^iogZie). (3.9) 

1/4 

Let us start with the macroscopic limit (1 <C V^hm? ~ V^hm?). In this limit the integral in 
(3.4) is dominated by contributions from saddle points. Depending on the sign of h+2h cos 0, 
the dominant saddle corresponds to either tiU = ±2 or trl7 = 0. With this taken into 
account, we obtain 

E(0) = -:! logZ^;i2(m,0) (3.10) 

~ —2m?{^gi + h + g 2 cos0 + \h + 2/icos0|) (3-11) 

up to subleading corrections. This function exhibits some interesting features. 

^^Recall that uo > 0 due to (3.6). 
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— B = 8 

•■■■ B = 6 

— B = 4 
■■■■ B = 2 
— B = 0 
— B = -2 

■■■■ B = -3 

— B = -4 

.... B = -6 

i.j z.u 

O/tt 

Figure 3. Dimensionless energy density (3.12) in the macroscopic limit for K = 4 and Nf = 2 
with equal masses at ^4 = g 2 /\h\ = 4 for varying B = 2h/h. 


/ When h > 2h, it follows (recall (3.6)) that h > 2\h\, so the energy becomes an 
analytic fnnction of O'. E{9) = —+ 2/i + {g 2 + 2h) cos0}. No phase transition 
is enconntered as 9 is varied. 

/ By contrast, when h < 2h, E{9) becomes non-analytic at those 9 where cos 9 = 
—h/2h. There are two first-order phase transitions in 0 < 0 < 27r. 

To examine the behavior of E{9) it is nsefnl to define the dimensionless energy density 

£{9) = — A cos 9 — \1 + B cos 9\ + 1 + A + B , (3.12) 

with A = <72/1^1 and B = 2hjh.^^ E{9) is plotted in fignre 3 for varying B ai A = A. We 
observe that E{9) has two cnsps for all \B\ > 1. At these first-order transition points, there 
are two degenerate vacna with eqnal energy density that can coexist by forming a domain 
wall. 

It is worth stressing that the 0-dependence of the energy density presented here dra¬ 
matically differs from that of QCD vacnnm. In the orthodox chiral perturbation theory, 
a first-order transition takes place at 0 = vr for two degenerate flavors and there CP is 
spontaneonsly broken [60, 61, 65, 67]. By contrast, nothing dramatic happens at 0 = tt in 
the Stern phase. 

An important remark on the topological charge distribntion is in order. When h > 
2\h\, ~ Q'iV 4 nA[gi+ 2 h+{g 2 + 2 h)cose] macroscopic limit. This means that 

topological charges are distribnted according to the weight 


Q \ / r^2\^ T t / r^2\^ exp [-QV(2 (Q^))] 

= -p[- (Q >] Iq((Q » ~-, 


(3.13) 


^^£{6) has the same 6 dependence as E{0)/2rrE\h\. The constant part of £{9) was chosen for an aesthetic 
reason (to ensure that curves of £{9) for different values of B in figure 3 do not intersect). 
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- 2V4|h|m2 = 0.1 

- 2V4|h|m- = 0.2 

- 2V4|h|m2 = 0.4 

- 2V4|h|m2 = 1 

- 2V4|li|m2 = 10 


Q jr: 


Figure 4. Finite-volume energy density (3.14) for K — 4 and Nf = 2 with equal masses for 
varying 2V4\h\m? at 2l = 4 and B = 8. In the limit 2V4\h\rnP' oo the curves converge to the 
B = 8 curve in figure 3. 


where = V 4 Xt° = 21 / 4771 ^( 5(2 + 2/i) and in the second step we have used an asymptotic 

formula for the modified Bessel function of first kind. Equation (3.13) is thus valid for 
1 <C (Q^) and Q <C (Q^)- Intriguingly, exactly the same functional form as (3.13) is known 
for the topological charge distribution in one-flavor QCD [46] and in high-temperature 
QCD [30, 68]; in both cases there are no massless Nambu-Goldstone modes because chiral 
symmetry is unbroken, and the topological charge obeys Poisson statistics. By contrast, 
the Stern phase do produce pions and yet exhibits the same topology dependence, which 
comes as a surprise. 

Next we leave the macroscopic limit and proceed to the finite-volume regime where 
microscopic variables take 0(1) values. This means that the zero-mode fluctuations of 
pions can no longer be ignored. The dimensionless finite-volume energy density can be 
defined, from (3.5), as 

£Rn{0) = -Acos9 - ^ — ^log{Io{ae) - hiae)} + 1 + A + B, (3.14) 

2 V 4 \h\m^ 

which reduces to £{0) in (3.12) as V4|/i|?77^ —>■ 00 . In figure 4 <£’fin(0) is plotted for various 
V 4 \h\m‘^ at fixed A and B. We observe that, while Tfin(^) is an analytic function of 6 , it 
gradually develops sharp peaks as V4|/i|r?7^ is increased. In the limit V 4 \h\m? —>■ 00 they 
turn into genuine first-order phase transitions, as depicted earlier in figure 3. 


• Topological density The 6 dependence of the vacuum can also be probed by the 
topological density ( ’ ’'^hich is defined in a dimensionless form as 


a{e) = 


d£{ 6 ) 

dB 


(3.15) 


in the macroscopic limit. This is plotted in figure 5. The discontinuous jumps of cj( 0) 
represent phase transitions. Intriguingly, for A = B = 4 there is a finite range of 9 where 
cr{9) vanishes exactly. This requires fine-tuning of low-energy constants and may not be 
realized in the real world, though. 
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Figure 5. Topological density in the macroscopic limit for if = 4 and Nj = 2 with equal masses 
at A = 4. 



B = 1 
B = 2 


- B = 4 

— B = 6 



0.0 0.5 1.0 1.5 2.0 


O/tt 

Figure 6. Phase diagram in the macroscopic limit for K = 4 and Nf = 2 with equal masses, 
with B = 2h/h. The energy is minimized by tr U = ±2 in the empty region and by tr ?7 = 0 in the 
shaded region, respectively. The phase transitions at the boundaries are generically first order. 

• Exotic flavor symmetry breaking The phase structure in the macroscopic limit is 
summarized in the phase diagram in figure 6, in the plane spanned by 9 and B = 2h/h. As 
can be seen from (3.4), phase transitions occur when h + 2hcos9 switches sign. Dividing it 
by h+2h (> 0; recall (3.6)) we get so the phase boundaries are set by 1+B cos 9 = 0 

and B = —1. In figure 6 we observe that the phase boundaries for |i?| S> 1 asymptote to 
9 = 7r/2 and 37r/2. This is because h + 2hcos9 ~ 2hcos9 changes sign at those 0’s. 

The phases in different colors in figure 6 exhibit distinctive properties. 

/ In the white region (h + 2hcos9 > 0), the energy is minimized at trI7 = ±2, i.e., 
U = ±1. One can parametrize fluctuations around 1 as 17 = exp (i(j)aTa/ f), insert 
this into (3.4) and expand the exponent up to second order in cjia, which enables us 
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S = 0.5 5 = 5 



Figure 7. Pion masses as a function of 9. In the left panel, all pions are degenerate. In the 
right panel, there are two branches that split at intermediate 6 . This is indicative of spontaneous 
breaking of SU(2)v. 


to read off the pion masses as 


m 


2 

TT 


Sm?{h + 2h cos 9) 

P 


The three pions are degenerate in this phase. 


(3.16) 


/ In the bine region (h + 2hcos9 < 0), the energy is minimized at tr t/ = 0 and leads 
to degenerate vacna. If we take U = fra as a representative and analyze qnadratic 
flnctnations aronnd it, the masses of three pions are fonnd to be 


m 


2 

TT 


0 , 0, and 


Sm?\h + 2hcos9\ 

P 


(3.17) 


The two gapless modes correspond to the vector rotations in 1- and 2-directions, while 
the gapped mode corresponds to the axial rotation in 3-direction. 


The mass spectra (3.16) and (3.17) are shown in hgnre 7 for B = 0.5 and 5. At i? = 5, 
two of the pions go massless at intermediate 0, signaling the breakdown of vectorial SU(2) 
symmetry. Since the vectorial flavor symmetry is an exact symmetry of QCD for degenerate 
masses, we expect that higher-order terms in the chiral effective theory wonld not spoil their 
masslessness, as long as the qnark masses are degenerate. 

Two remarks on the literatnre are in order. First, the exotic flavor-breaking phase fonnd 
above is analogons to the so-called Aoki phase [17] in lattice QCD with Wilson fermions. It 
has been shown by Sharpe and Singleton that the Aoki phase originates from a competition 
among terms at 0{m), 0(a) and 0(P) in the chiral Lagrangian [70], with a the lattice 
spacing. This is similar in essence to onr effective theory for the Stern phase, in which 
competing terms arise at 0{rP). Secondly, it has been pointed ont by Crentz [36] and 
Smilga [37] for Nf = 2 and 0 ~ tt that a similar vectorial flavor breaking can take place 
even in the standard chiral effective theory if a particnlar sign is chosen for a low-energy 

^^Note that the Vafa-Witten theorem [69] prohibiting vectorial symmetry breaking does not apply at 
finite baryon density or nonzero 0 angle. 
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constant at 0{p‘^). Although their analysis has nothing to do with the Stern phase, the 
technical aspects of their analysis are similar to ours. 

The extension of results in this section to non-degenerate masses or to > 2 would 
be technically more involved. This is deferred to future work. 


3.2 Quarks in higher representations 

In this section we consider 0-dependence of QCD-like theories with Nf > 1 flavors of Dirac 
fermions in a general complex representation R of the gauge group. The motivation for such 
an extension comes from several directions. First, gauge theories with fermions in higher 
representations have attracted interests as promising candidates of the beyond-Standard- 
Model physics [71-73]. Secondly, large-QCD with quarks in the adjoint, two-index 
symmetric and antisymmetric representations of the gauge group are of interest from the 
viewpoint of orientifold planar equivalence [74, 75]. Of course, whether the Stern phase can 
be realized in such theories is a highly nontrivial dynamical question for which we have no 
definitive answer yet. In what follows, we shall take the existence of the Stern phase as an 
assumption and discuss outcomes specific to quarks in higher representation. 

First and foremost, the index theorem states that Iji = 2 TrQ in the background of 
gauge fields with the topological charge Q. Then the 0-angle enters the partition function 
only through the combination e*®(det or in other words, the 0 dependence can be 

incorporated into effective theory via an axial rotation 

Looking back at (2.14), we again find that the case with K > A has no 0 dependence at 
leading order of the £ expansion. To see 0-dependent physics and for the sake of technical 
simplicity, we concentrate on the Nf = 2 and K = A case in the following. Furthermore, 
to make the discussion explicit, we will take R to be the sextet (two-index symmetric) 
representation of SU(3),^® for which Tr = 5/2, although any other higher representation 
will do the job. From (2.4) the non-anomalous subgroup of U(l)y 4 is Z 20 , which is supposed 
to be spontaneously broken to J-k = ^ 4 - 

Plugging M = into (2.14b) one finds the finite-volume partition function for 

sextet fermions. 


z§J2{m,e) = 


dU exp 


/SU(2) 

+ (^2e*®/5+h.c.)} 


V 4 m^|/i(tr t/)^ -|- 2 gi + (he*®/® tr(t/^) -|- h.c.) 


(3.19) 


A new interesting feature of this partition function is that it is periodic in 0 with period 
lOvr, rather than 27r. This appears to contradict the 27r-periodicity of (3.1). The resolution 
of this “puzzle” goes as follows. As noted above, the theory with sextet quarks has Z 20 
unbroken axial symmetry in the chiral limit. The putative higher-order (e.g., quartic) 
quark condensate is invariant only under Z 4 C Z 20 , so there are five degenerate vacua. (We 

*®Asymptotic freedom requires Nf < 3. 
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remind the reader that the existence of hve isolated components of the vacuum manifold 
follows from (2.6) in section 2.1.1.) Once we switch on the quark masses, the hve-fold 
degeneracy is lifted and one of those vacua is selected as the unique ground state. In fact, 
the effective theory (3.19) is a theory of fluctuations around such a ground state. Now, if 
we rotate the 9 angle gradually, those hve vacua are permutated in a cyclic way and the 
ground state moves from one state to another. After a 27r rotation of 6 , those hve low-lying 
states undergo a cyclic rotation by one unit, and the system as a whole returns to itself, 
despite that each state returns to itself only after lOvr rotation of 9. 

We now have two comments: 

• This mechanism was already pointed out by Leutwyler and Smilga [46] for SU(A^) 
gauge theory with adjoint quarks. They explained how the 27rAi-periodicity of the 
effective theory in 9 can be reconciled with the 27r-periodicity of the full theory. As 
noted in [76], this also pertains to the well-known subtlety that the pure Yang-Mills 
partition function can be 27r-periodic in 9 even though the large-Ai scaling tells that 
the natural variable in the large-limit is 9/N rather than 9 [60, 77]. 

• When quarks are in the fundamental representation (T^ = 1/2), the vacuum man¬ 
ifold in the chiral limit only has a single connected component [cf. (2.6) in section 
2.1.1]. This means that in the case of fundamental quarks we need not sum up contri¬ 
butions from multiple disconnected sectors explicitly to recover 27r-periodicity of the 
full partition function. 

The full partition function for the Stern phase with sextet quarks may be dehned as 

1 ^ 

^m(0) - 5 E 9 + 27rk) , (3.20) 

k=0 

which is manifestly 27r-periodic in 0. If we take the macroscopic limit in the e-regime, the 
state having the lowest energy will dominate (3.20). The energy density is therefore 

E^{9) := min E (3-21) 

0<A:<4 \ 5 J 

with E{9) in (3.11). The dimensionless energy density can be dehned similarly as 

^^(0) := min £ {3.22) 

o<fc<4 \ 5 J 

with £{9) in (3.12). Figure 8 displays £’m(^) for 0 < 9 < lOvr with two sets of A and B. One 
can clearly observe the branched structure of the 9 vacuum. In the left panel, hrst-order 
phase transitions occur at 0 = {2i + l)7r with ^ ~L. Notably, the exotic havor-breaking 
phase found for fundamental quarks in section 3.1 does not appear. This is because for 
A = A and B = 8 the exotic phase tends to have higher energy than the normal phase and 
is consequently disfavored in the minimization in (3.22). 

However this is not necessarily true for other values of A and B. By decreasing A one 
can lower the energy of the exotic phase at will. As an example, we show in the right panel 
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Figure 8. Dimensionless energy density (3.22) [thick solid line in orange] in the macroscopic limit 
for iir = 4 and Nf = 2 with equal masses in the sextet representation. The k = 0 branch is denoted 
by a thick blue dashed line. 


of figure 8 the energy density for A = —2 and 5 = 3. In this case the first-order phase 
transitions occur at 6 = 2 i 7 r with I At 0 = 0 the vacuum is two-fold degenerate and 
breaks parity spontaneously.^® Moreover, we discover that the exotic flavor-breaking phase 
is realized for all values of 6 This tells us that the condition (3.6) is not sufficient, in the 
case of quarks in higher representations, to ensure that U = 1 is the ground state at 0 = 0. 
In short, the Stern phase with quarks in higher representations exhibits an “all-or-nothing” 
behavior: if the vectorial flavor symmetry is unbroken at 0 = 0 , it is unbroken at any 9, 
and conversely, if it is broken at 0 = 0, it remains broken at any 9. Since this is possible 
only if multiple states exchange dominance for varying 9, it cannot happen for quarks in 
the fundamental representation that have only one vacuum sector. 

4 Conclusion 

We have investigated properties of the Stern phase using the low-energy effective theory of 
pions at zero and nonzero vacuum angle 9. Analytical results are obtained for the 9 and 
volume-dependence of miscellaneous physical quantities, both for fundamental quarks and 
for quarks in higher representations of the gauge group. We have highlighted an intricate 

^^Spontaneous parity breaking at finite chemical potential is not ruled out by the Vafa-Witten theorem 
[78] because the path-integral measure becomes complex. We refer to [79-81] for recent model studies 
of parity-breaking phases at finite density. It also deserves attention that the Vafa-Witten theorem [78] 
can fail for quark bilinears even with positive-definite measures [17, 82] ; various authors have investigated 
limitations of the original proof of the theorem [70, 83-87]. It would be quite interesting to extend the 
proof so as to incorporate four-fermion condensates, though we do not attempt it here. 

^^This can be checked as follows. Let us first notice that the fc = 0 branch is the ground state for 
4% < d < Ott. This implies that, for 6 in this range, it is the sign ot h 2/icos ( ^+g’”'*' ) = h 2hcos(0/5) 
that determines whether the exotic flavor-breaking phase is realized or not. Since B + 1 = {h2h)/h = 
4 > 0 we have h > 0. Next note that 1 + Bcos{9/5) < 0 for 9 £ [47r,67r]. These together imply that 
h 2/icos(0/5) = /i{l -f Bcos{9/5)} < 0. Therefore the flavor-breaking phase with trU = 0 is realized. 
Upon inspection this is seen to extend to all 9. 
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interplay of multiple competing terms in the chiral Lagrangian and discussed its relevance 


for the phase structure of the Stern phase at nonzero 6 . Instead of the standard Dashen’s 


phenomenon at 0 = tt, we have found either two first-order phase transitions or no transition 
at all, depending on the values of low-energy constants. Throughout this work we have only 
relied on symmetries of the system. Therefore the obtained results should be robust as long 
as a nonzero mass gap exists for non-Nambu-Goldstone modes. 

Since lattice simulations at finite density or 0 7 ^ 0 are currently unfeasible, it will be 
worthwhile to extend and improve theoretical examinations along the lines of this work 
further. There are several future directions. Firstly, we can generalize our analysis of the 
9 vacua in section 3 from Nf = 2 to more flavors. Since there are more independent terms 
at 0{M‘^) than for Nf = 2 we can expect richer physics. A thorough study of profiles 
of domain walls discussed in section 2.1.1 may also be intriguing. Throughout this work, 
we did not attempt to find out the microscopic mechanism that realizes the Stern phase 
in finite-density QCD. This is a challenging open problem that no doubt deserves further 
investigation. Another interesting direction is to extend the present work to QCD-like 
theories with quarks in (pseudo)real representations of the gauge group. These theories 
enjoy extended flavor symmetries and it is interesting to ask how to define the Stern phase 
in this case. From a phenomenological point of view it is important to incorporate the effects 
of isospin chemical potential into the effective theory, which can be done along the lines of 
[ 88 , 89]. Analytical calculation of the unitary integrals in (2.14) is an open mathematical 
problem. 
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A Derivation of Xud in (2.24) 

This appendix outlines the derivation of (2.24). For brevity we introduce a shorthand 
notation for a group average: 



SU(2) 


dU f{U ). 


(AT) 
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Then 


dniudm*^ 


= Vih- 


= Vih 


'^\og(^(^eV4htr(MU)tr(WMr)\^\^ 

Q {(tx{MU) U^2 

^^e'^4/itr(Mf/)tr(t/tMt)^^ 


U\\U* 0^4^ tr(MC/) tr(C/tMt) 


+ {Vihf 


- {V^hf 


^^gy4/itr(M!7)tr((7tMt)^^ '' 

IUt{MU) [/|2e^4/itr(Mt/)tr(C/tMt) 


UiiU^2 tT{MU) tr(C/tMt) eV4/itr(Mf/)tr(t/tMt) 
^^0^4/1 tr(M! 7 )tr(( 7 tMt)^^ 

('tr([/tMt) C/n e^4/itr(MC/)tr(C/tMt)\\ 


^^0^4/1 tr(M( 7 ) tr(f/tMt) ^ 


After switching to the microscopic variables Hf = 2\/V^ nif and taking the degenerate 
mass limit, we get 

n P 4 omuOm^ \\ II 




,i^(trC/)2 


log 




^2 ^^(trl7)C/2*2e'^^^''^^"^^ ^(|(trC/)C/iie^(*''^)' 


To compute the group average it is convenient to adopt the parametrization based on 
SU(2) ^ 5^: 


Xo + iX 3 X 2 + ixi 
-X 2 + ixi Xo - iX 3 


with Xg + X 4 + X 2 + x| = 1 , 


Then the followings hold for arbitrary a G IR: 

((02-g^^ =0«{/g(a)-/i(a)} , (A.6a) 

((x^2“"o^) = l0“{/g(a) -/2(a)} , (A.6b) 

((xie2“"o^) = le“{3/o(a)-4/i(a) + /2(a)} , (A.6c) 

((x^e2“"o^) = le“{2/o(a) + h{a) - 2l2{a) - h{a)} , (A.6d) 

((xgx|e2""^o)) = ^e"{2/o(a) - Ii(a) - 2/2(0) +/3(a)} , (A.6e) 

where In{(^) is the modihed Bessel function of the hrst kind. 
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We now substitute (A.4) into (2.23) and use (A.6) with a = which yields 


Xud 





I h-h + ^ ( lo-h V 

^ lo — h 12 Iq — h 16 V/q —1^1/ ’ 







(A.7) 

(A.8) 


with the argument /i^/2 omitted. This is the desired result. 
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